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ABSTRACT 



Aims. Off-the-limb observations with high spatial and spectral resolution will help us understand the physical properties of spicules 
in the solar chromosphere. 

Methods. Spectropolarimetric observations of spicules in the Hei 10830 A multiplet were obtained with the Tenerife Infrared 
Polarimeter on the German Vacuum Tower Telescope at the Observatorio del Teide (Tenerife, Spain). The analysis shows the variation 
of the off-limb emission profiles as a function of the distance to the visible solar limb. The ratio between the intensities of the blue 
and the red components of this triplet (ft = /biue/^red) is an observational signature of the optical thickness along the light path, which 
is related to the intensity of the coronal irradiation. 

Results. We present observations of the intensity profiles of spicules above a quiet Sun region. The observable ft as a function of 
the distance to the visible limb is also given. We have compared our observational results to the intensity ratio obtained from de- 
tailed radiative transfer calculations in semi-empirical models of the solar atmosphere assuming spherical geometry. The agreement is 
purely qualitative. We argue that future models of the solar chromosphere and transition region should account for the observational 
constraints presented here. 
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1. Introduction 

The solar chromosphere is in a highly dynamic state. Its small- 
scale structures evolve in timescales of minutes or even less. Due 
to its low density, the chromosphere is transparent in most of 
the optical spectrum. Nevertheless, in lines such as Ha and He 
i 10830 A, the significant absorption provides a mean for direct 
studies of the chromosphere's peculiar characteristics, such as 
bright plages in active regions, dark and bright fibrils on the disc 
in both the quiet Sun's network and around su nspots, as well 
as spi cules above the limb. Recent studies (e.g. iTziotziou et alj 
120031) claimed that it is possible that many of these chromo- 
spheric features have the same physical nature within different 
scenarios. 

In the chromosphere, between the photosphere and the much 
hotter corona, the average temperature starts to rise outwards, 
along with a decrease of the density. The dynamics of a mag- 
netised gas depends on the ratio of the gas pressure and mag- 
netic pressure, the plasma B parameter, B-P g:ls /P mag , with 
Pmag = B 2 /(87r). From the low chromosphere into the extended 
corona, this plasma parameter decreases from values B > 1 to a 
low -beta regime, B <k 1, where the plasma motions are magnet- 
ically driven. 

Spi cules, known for more than 130 years (see hand draw- 
ings in ISecchill 18771) . represent a prominent example of the dy- 
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namic chro mosphere. We refer th e reader to review s bv lBeckersI 
dl968lll972l) and to the paper bv lWilhelml d2000l) on UV prop- 
erties. Spicules are seen at and outside the limb of the Sun as 
thin, elongated features that develop speeds of 10-30 km s _1 and 
reach heights of 5-9 Mm on av erage, during th eir lifetimes of 
3-15 minutes. As pointed out bv lSterli ng (2000), a key impedi- 
ment to develop a satisfactory understanding has been the lack of 
reliable observational data. Many theoretical models have been 
developed to understand the nature of spicules, using a wide va- 
riety of motion triggers and driving mechanisms. In this study 
we focus on the Hei 10830 A triplet emission line using recent 
technical improvements in observational facilities. We are able 
to provide observational evidence of the link between the corona 
and the infrared emission of this line, in the frame of the current 
theoretical models of the solar atmosphere. 

The Hei 10830 A multiplet consists of the three transitions 
from the upper term 3 P2,i,o, which has three levels, to the lower 
metastable term 3 Si, which has one single level. The two transi- 
tions from the J=2 and J= 1 upper levels appear blended at typical 
chromospheric temperatures, and form the so-called red compo- 
nent, at 10830.3 A. (Note that the two red transitions are only 
0.09 A apart). The blue component, at 10829.1 A, corresponds 
to the transition from the upper level with J=0 to the lower level 
with J=l. The energy levels that take part in these transitions 
are basically popul ated via an ionization-recombination process 
dAvrettet alJll994l) . The EUV coronal irradiation (CI) at wave- 
lengths lambda A < 504 A ionizes the neutral helium, and sub- 
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sequent recombinations of singly ionized helium with free elec- 
trons lead to an overpopulation of these levels. 

Alternative theories suggest other mechanisms that might 
also contribute to the formation of the helium lines relying on 
the collisional excit ation of the electrons in regions with higher 
temperature (e.g jAndretta & Jonesll 1997b 

ICentenol d2006h modelled the ionisation and recombination 
processes using various amounts of CI, non-LTE radiative trans- 
fer, a nd different atmospheric models (see also ICenteno et al] 
120071 in preparation). They have simulated limb observations for 
different heights, obtaining synthetic emission profiles in spheri- 
cally symmetric models of the solar atmosphere. One conclusion 
of their study is that the ratio of intensities CR = /biue/^red) of the 
'blue' to the 'red' components of the He 1 10830 A emission is a 
very good candidate for diagnosing the CI. The population of the 
metaestable level depends on optical thickness, whose variation 
with height governs the change in the ratio 'R as a function of the 
distance to the limb. 

Truiil lo Bueno et ail d2005l) measured the four Stokes pa- 
rameters of quiet-Sun chromospheric spicules and could show 
evidence of the Hanle effect by the action of inclined magnetic 
fields with an average strength of the order of 10 G. They mod- 
elled the He i 10830 A profiles assuming the medium along the 
integrated line of sight as a slab of constant properties and with 
its op tical thickness as a free parameter. Trui illo Bueno et aT] 
(2005) showed that the observed intensity profiles and their 
ensuing 'R values can be reproduced by choosing an optical 
thickness significantly larger than unity. ICentenol d2006l) demon- 
strated that this optical thickness is related to the coronal irra- 
diance (through the ratio K), thus providing a physical meaning 
to the free parameter in the slab model (see also Centeno et al. 
2007). 

We present novel observations showing the spectral emission 
of He i 10830 A and its dependence on the height of the spicules 
above a quiet region. We compare the deduced observational 'R 
with that obtained from detailed non-LTE numerical calculations 
using available atmospheric profiles. 



2. Observations 

The observations were carried out on December 4 rh , 2005, at the 
Vacuum Tower Telescope (VTT) at Observatorio del Teide. They 
were sup ported by the Kiepenhe uer Adaptive Optics System 
(KAOS. [von der Liihe et alj |2003). We used the Echelle spec- 
trograph of the VTT a nd the new version (TIP-2) of the Tenerife 
Infrared Polarimeter (Martin ez Pillet et"al]|1999t) . which has a 
larger CCD camera dCollados et alj 120071) . The seeing condi- 
tions were good (average seeing after KAOS correction around 
7 cm, maximum 12 cm). 

The strong darkening close to the solar limb and the actual 
presence of the limb make it difficult to use KAOS for off-limb 
observations, since the correlation algorithm was not developed 
for this kind of observation. Fortunately, we could use a nearby 
facula inside the limb for the tracking lock point of the adaptive 
optics system. 

Although we carried out full Stokes vector observations, in 
this study we have used only the intensity profiles. The spectral 
region covered by TIP-2 spans from 10826 to 10837 A and con- 
tains the He i 10830 A multiplet. The spectral sampling was 11 
mA/px, and the slit was 40" long and 0'.'5 wide. It was oriented 
parallel to the NE limb, in a region far from major activity at 
that time, at 59?8 from North. We scanned the full height of the 



spicules, up to 7" off the visible solar limb, with a step size of 
a/35. 

At each position, 5 consecutive spectra were measured with 
10 accumulations of 250 ms each, with a total integration time 
of 2.5 s per off-limb slit position. A nearby disc profile was also 
taken to help remove the scattered light from the off-limb spectra 
via the data analysis. 



3. Data reduction and processing 

We applied usual flatfielding and dark current corrections to the 
data. A high frequency, quasi-periodic, spurious electronic pat- 
tern in the profiles was removed using a low-pass Fourier filter, 
which left the frequencies containing the spectral line informa- 
tion untouched. 

We define the position of the solar limb as the height of the 
first scanning position where the helium line appears in emis- 
sion. For increasing distances to the solar limb a decreasing 
amount of sunlight is added by scattering in the Earth's atmo- 
sphere and by the telescope optical surfaces to the signal. Since 
the true off-limb continuum must be close to zero, i.e. below our 
detection limit, the observed continuum signal measures the spu- 
rious light. Therefore, we removed the spurious continuum in- 
tensity level by using the information given by a nearby average 
disc spectrogram. This first subtraction estimates the continuum 
level on a region 6 A away from the Hei 10830 A emission 
lines. After this correction with a coarse estimate of the spuri- 
ous light, a second correction was applied to remove the residual 
continuum level seen around the emission lines. This was needed 
since the transmission curve of the used prefilter is not flat but 
variable with wavelength. 

FigureQ]shows the emission profiles of the He 1 10830 A (af- 
ter the reduction process) for different heights above the limb. 
Figure [2] illustrates this in three dimensions, as a function of 
wavelength and the distance to the solar limb, clearly showing 
a change in the intensity ratio of the blue and red components of 
the multiplet CR = /biue/^-ed) with height. This will be discussed 
in Sect. [4] 
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Fig. 1. Measured He i 10830 A emission profiles for increasing 
distances to the solar limb, scanning a broad range of the height 
extension of the spicules. Each profile is the average of the 312 
pixels along the slit (which was always kept parallel to the limb) 
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Fig. 2. 3D representation of the measured He 1 10830 A emission 
profiles for increasing distances to the solar visible limb. Note 
that the x-axis is wavelength, the y-axis the height above solar 
limb and the z-axis the intensity normalised to the maximum 
emission in the line centre of the red component. 



For the calculation of % we need to determine the amplitudes 
of the blue and red components of the emission profile (as shown 
in Fig.[3]l. 

To determine the core wavelength of the red component of 
the triplet we fitted a Gaussian profile to its core, in a 1.3 A 
range around the maximum. After symmetrising the observed 
profile around this maximum, using the values on the red side 
of the red component, we fitted another Gaussian function to 
the resulting symmetric profile. Subtraction of the fitted sym- 
metric profile from the data leaves the emission profile of the 
blue component, which was also approximated by a Gaussian to 
determine its central wavelength. Our tests trying to fit directly 
both profiles using two Gaussians failed in a number of cases, 
probably due to the following reasons: (a) the red component is 
in fact the result of two blended lines, (b) the much weaker blue 
component was almost hidden in the broadened red component, 
and (c) the presence of noise. Our technique determines first the 
red component and then, after subtraction of the fitted profile, 
the blue one. 

We have thus separated the helium emissions into their red 
and blue components assuming only 2 that both are present and 
that they are both symmetric. We can now measure their widths 
and intensities and also check that the line core positions coin- 
cide with the theoretical ones. After the fitting process the resid- 
uals between measured and observed profiles were small, the 
largest errors occurring in the determination of the core inten- 
sities of the red line. This happens because the red component 
consists of two blended lines (with a separation of 0.09 A), a 
fact that flattens the emission profile near the core as opposed 
to a more peaked Gaussian function. Nevertheless, the differ- 
ences between fitted profiles and data are only significant in the 
red core and are always lower than +0.08 of the maximum nor- 
malised intensity, with a mean difference of ~0.02. To avoid sys- 
tematic errors, we used the observational values for the centre of 
the red component when calculating ft. 

4. Results 

The chromospheric temperature and density are too low to pop - 
ulate the ortho-helium levels via collisions dAvrett et al.| [l994). 
The EUV irradiation from the corona (CI) ionises the para- 
helium, and the subsequent recombinations lead to an overpopu- 
lation of all the ortho-helium levels, in particular those involved 




10828. 10829. 10830. 10831. 10832. 10833. 

A 

Fig. 3. Determination of the blue and red components of the He i 
10830 A triplet from the observed emission profiles. In this ex- 
ample the slit was placed at 1'.'4 off the solar visible limb. See 
text for details. The solid line represents the average emision 
profile. The dotted line is the Gaussian fit to the symmetrised red 
component. Subtraction of this from the observed profile leaves 
the blue component, which is also fitted by a Gaussian profile 
(thin solid line). The sum of both Gaussians (dashed line) gives 
the fit to the observed profile. 



in the 10830 A transitions. ICentenol d2006) and Centeno et al. 
(2007) have modelled the off-the-limb emission profiles and 
concluded that the ratio H = I blue 1 1 red is a function of the height 
and a direct tracer of the amount of CI. Here we compare the 
results from the theoretical modelling with observations. 

Trujill o Bueno et alj d2005l) modelled their spectropolari- 
metric observations assuming a slab with constant physical prop- 
erties with a given optical thickness. In the optically thin regime 
!R~0.12, which is the ratio of the relative oscillator strengths 
of the triplet. As the optical thickness (at the line-centre of the 
red blended component) grows, this ratio also increases until 
it reaches a saturation value slightly larger than 1 for r ~ 10. 
(This type of calculation can be done and improved as explained 
in Trujillo Bueno & Asensio Ramos 2007). To repr oduce the 
observed emission profile Truji llo Bueno et all (2005) had to 
choose t ~ 3. Interestingly, the values of t yielded by this mod- 
elling strategy are consistent with the more realistic approach of 
Centeno (2006), where non-LTE radiative transfer calculations 
in semi-empirical models of the solar atmosphere are presented, 
using spherical geometry and taking into account the ionising 
coronal irradiation. With our data we are able to test such the- 
oretical calculations by comparing the measured values of % 
with those resulting from various chromospheric models. This 
way we may eventually trace the amount of CI inciding on the 
spicules. The analysis described in Sect. [3] yielded the values 
of H for the observed profiles. The resulting dependence on the 
distance to the solar limb, for each pixel along the slit and each 
position of the slit above the limb, are presented in Fig. [4] The 
solid black line gives the average value of fl. 

The dependence of % with height can be understood in a 
qualitative way as follows: in the outer layers of the chromo- 
sphere the density is so low that the transitions occur in the opti- 
cally thin regime. With decreasing altitude the ratio % increases 
(proportionally with density) until a maximum optical thickness 
is reached. At even lower layers, although the density still con- 
tinues to rise, the extinction of the coronal irradiance leads to 
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Fig. 4. Measured ratio K- 1 blue 1 1 red as a function of distance to 
the solar limb. Thin lines come from each pixel along the slit. 
The thick solid line represents the average and the dashed line 
the value of the optically thin regime. 



a reduction in the number of ionizations, which results in a de- 
crease of the optical thickness in the core wavelength of the red 
component, and thus in a decrease of %. 

For a quantitative compari son with theor etical modelling we 
have u sed the results from ICentenol y006) and ICenteno et alj 
(2007) where they calculated the ratios H for different standard 
model atmo spheres: FAL -C and FAL-P dFontenla et al J 1 1 99 lb 
and FAL-X dAvretjl 19951) . The FAL-C and FAL-X models may 
be considered as illustrative of the thermal conditions in the 
quiet Sun, while the FAL-P model of a plage region. The FAL-X 
model has a relatively cool atmosphere in order to explain the 
molecular CO absorption at 4.6 jj.. 

The comparison is shown in Fig. [5] We notice that the mod- 
elled height variations of % agree only in a qualitative manner 
with what is found in our observations. However, the calcula- 
tions from different models of the solar atmosphere are unable 
to reproduce the measured ratio. Higher values of the coronal ir- 
radiance lead to an increase of the optical thickness (at the line 
centres of the He i multiplet) and an upward shift in the run of 
% vs. height. Yet the shape of the height dependence is mainly 
given by the atmospheric density profile and the attenuation of 
the ionising radiation as it reaches the lower layers of the chro- 
mosphere. It is also clear from Fig. [5] that the models do not 
extend high enough. 



5. Conclusions 

The theoretical behaviour of the ratio % agrees qualitatively with 
observations. Yet, a quantitative comparison shows poor agree- 
ment. Also, the simulated ratios are highly model dependent. As 
already explained, the failure to reproduce the observed pro- 
files is very likely due to the density stratification not being ade- 
quate for spicule modelling and to the limited vertical extension 
of the atmospheric models. Modelling of the intensity ratio % in 
the He i infrared triplet should account for the fact that the so- 
lar chromosphere is inhomogeneous on small scales and that the 
spicules are small-scale intrusions of chromospheric matter into 
the hot corona. 

New data of spicule regions near the poles and the equator, 
below coronal holes or coronal active regions should help us to 
understand the detailed behaviour of the Hei 10830 A lines. In 
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further work, we will extend this study to the full Stokes vector, 
in order to see the variation of the linear polarization - or even 
the variation of the Hanle effect - with height. It would also be 
interesting to use the most recent mode ls of active region fibrils 
and spicules (e.g jHeggland et ai1l2007l) in order to see whether 
or not they agree with our observations. Future models of the 
solar chromosphere should be constrained by the observational 
evidences presented here. 
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